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ABSTRACT 

Context. Spectroscopic observations and theoretical models suggest resonant wave-particle interactions, involving high-frequency 
ion-cyclotron waves, as the principal mechanism for heating and accelerating ions in the open coronal holes. However, the 
mechanism responsible for the generation of the ion-cyclotron waves remains unclear. One possible scenario is that ion beams 
originating from small-scale reconnection events can drive micro-instabilities that constitute a possible source for the excitation 
of ion-cyclotron waves. 

Aims. In order to study ion beam-driven electromagnetic instabilities, the multi-fluid model in the low-/3 coronal plasma is 
used. While neglecting the electron inertia this model allows one to take into account ion-cyclotron wave effects that are absent 
from the one-fluid MHD model. Realistic models of density and temperature as well as a 2-D analytical magnetic field model 
are used to define the background plasma in the open-field funnel region of a polar coronal hole. 

Methods. Considering the WKB approximation, a Fourier plane-wave linear mode analysis is employed in order to derive the 
dispersion relation. Ray-tracing theory is used to compute the ray path of the unstable wave as well as the evolution of the 
growth rate of the wave while propagating in the coronal funnel. 

Results. We demonstrate that, in typical coronal holes conditions and assuming realistic values of the beam velocity, the 
free energy provided by the ion beam propagating parallel the ambient field can drive micro-instabilities through resonant 
ion-cyclotron excitation. 
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1. Introduction 

Observations made by the Ultraviolet Coronagraph 
Spectrometer (UVCS) and other instruments on the Solar 
and Heliospheric Observatory (SOHO) have significantly 
increased our knowledge of the kinetic properties of 
charged particles close to the Sun in the source region 
of the fast solar wind. Spectroscopic determination of the 
widths of ultraviolet emission lines in coronal holes indi- 
cate that heavy ions are very hot and have high temper- 
ature anisotropies, and that heavier ions have a higher 
temperature than the protons by at least their mass ra- 
tio, i.e. Ti/Tp > mi/nip ( Kohl et al.||1997 Cranmer et al 



1999[). These observations strongly suggest resonant ion- 



cyclotron wave-particle interaction as a major mechanism 
for the heating and acceleration of ions in the magnetically 
open corona. This notion led to a renewed interest in mod- 
els involving ion heating by high-frequency ion-cyclotron 



waves (Isenberg et al. 2000 Hollweg 2000 Marsch fc Tu 



2001 Vocks & Marsch 2001 Xie et al. 2004). For a de- 



tailed review on resonant ion-cyclotron interactions in the 
corona, see 



Hollweg fc Isenberg (2002) 



However, it remains unclear how these waves originate 
in the solar corona, and whether they are generated locally 
or emanate from the coronal base. One possible scenario 
is that the ion-cyclotron waves are generated in the lower 
corona by small-scale reconnection events as suggested by 



Axford & McKenzie (1992). They might also be gener- 
ated locally through a turbulent cascade of low-frequency 
MHD-type waves towards high-frequency ion-cyclotron 



waves ( 


Li et al. 


1999 


Hollweg 


2000 



or by plasma microinstabilities driven by current fluctua- 
tions of low- frequency MHD modes ( Markovskii|2001 1 , or 
by an intermittent electron heat flux accompanying mi- 
croflare events (Markovskii & Hollweg 2004). Voitenko & 
Goossens (20021 further suggested that plasma outflows 
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from reconnection sites in microflares could create ion- 
beam configurations in the surrounding plasma, and thus 
provide free energy for driving kinetic microinstabilities 
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through ion-cyclotron resonance and the Cerenkov effect. 
The possible origin of the ion beams observed in the solar 
wind from reconnection jets and explosive events in the 



corona has been proposed by Feldman et al. (19961 



Indeed, detailed spectroscopic studies of the so-called 
high- velocity events and explosive events, using spec- 
tra obtained with the Coronal Diagnostic Spectrometer 
(CDS; |Brekke et ah] ( |1997| )) and the Solar Ultraviolet 
Measurement of Emitted Radiation instrument (SUMER; 
Innes et al. ( 1997[ )) both on SOHO, revealed a new charac- 
ter of the lower corona as a highly dynamic medium. They 
signify the omni-presence of transient explosive events and 
a wide variety of plasma jets with velocities ranging from a 
few tens of a kilometer per second up to several hundreds 
of kilometers per second. Since these plasma jets have been 
observed to evolve in a similar way as predicted by the 




Fig. 1. Magnetic field geometry of a funnel as obtained 



from the 2-D potential-field model derived by Hackenberg 



theory of magnetic reconnection (Innes et al. 1997), ex- 



plosive events and plasma jets have been associated with 
the highly-dynamic small-scale reconnections which are 
supposed to take place in the chromospheric network, ap- 
proximately at heights of 1000 — 3000 km above the pho- 
tosphere. Solar magnetograms provided by the Michclson 
Dopplcr Imager (MDI) on SOHO clearly indicate that the 
magnetic network field exists in two characteristic com- 
ponents, i.e. in few unipolar open flux tubes, "funnels", 



et al. (2000). The field lines emerge from the boundary 
between two adjacent supergranules (x — 0) and expand 
rapidly to fill the corona. The photospheric level is at z — 
0. Coordinate axes, wave vector and beam drift velocity 
are shown in red color. 



cussed in Sect. |4j and finally we give our conclusions in 
Sect. 02 



and in multiple closed flux tubes, "loops". (Schrijver et al. 2. Background plasma configuration 



19981. The magnetic network is very dynamic and re- 



leases non-potential-field magnetic energy which is con- 
verted into plasma heating, beam particles, and motion of 
the coronal plasma. All these processes provide ample free 
energy for driving plasma macro- and micro-instabilities. 

The present paper aims at studying the possible occur- 
rence of these instabilities in typical conditions for coronal 
holes. We particularly focus our study on beam-driven in- 
stabilities in the ion-cyclotron frequency range which are 
driven by the presence of tenuous ion beams, presumably 
originating from small-scale reconnection sites and propa- 
gating parallel to the ambient magnetic field. Linear mode 
analysis is used in the framework of the multi-fluid model, 
while neglecting the electron inertia. This model permits 
the consideration of ion-cyclotron-wave effects that are 
absent from the one-fluid MHD model. Realistic models 
of the density and temperature, as well as a 2-D funnel 
model describing the open-field region, are used to define 
the background plasma. Considering the WKB approxi- 
mation (in which the wavelength of interest is assumed 
to be smaller than the non-uniformity length scale), we 
first solve the dispersion relation locally and then perform 
a non-local wave analysis using the ray-tracing theory, 
which allows us to compute the ray paths of the unsta- 
ble waves in the funnel as well as the spatial variation of 
their growth rates. 

This paper is structured as follows. In Sect. [2] we 
present the 2-D analytical funnel model used in this study 
to describe open-field region in a coronal hole. Then in 
Sect. [3] we describe how the local and non-local (ray 
tracing) linear perturbation analysis using the multi-fluid 
model is carried out. The results are presented and dis- 



For the background plasma density and temperature we 



use the model parameters of Fontcnla et al. (1993) for the 



chromosphere and Gabriel ( 1976 ) for the lower corona 



The 2-D potential-field model derived by Hackenberg et al.| 
(2000) is used here to define the background magnetic field 
(Fig. [ll in the funnel. Analytically, the two components 
of this model field are given by: 



{B max ~B m )L 1 cosh^ ~ cos(f + 
B 0x (x, z) = — ~ 7-^ « In- 



2tt(L - d) cosh 



2txz 



cos{ 



Bq z (x, z)—B m + (B rnax — B t 



no) 



7rd 
L 



L 



2lTX \ 

L I 



(1) 



arctan- 



^ 2ttz ol ;^ nd 



L-d {L-d)ix 

L ) , 



sm(f£ 1 27rx ' 



arctan- 



cosh- 



sm 



2L 



, 27TZ 




(2) 



The typical parameters relevant for this model are: L — 
30 Mm, d = 0.34 Mm, B 00 = 11.8 G, and B max = 1.5 kG. 

3. Linear perturbation analysis 

To describe wave propagation in the funnel we use the 
multi-fluid equations and subject them to a linear pertur- 
bation analysis. The fluid equations associated with the 
polytropic gas law for any particle species j are given by: 



drij 
dt 



+ V ■ {rijVj) = 0, 



(3) 
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mjrij(-^- + Vj ■ Vvj) + Vpj — qjrij(E + Vj x B) 



X! " !v -' v i') = °> 



(4) 



Pjn- 7j = const, 



where nij, n 



., ( 5 ) 
j , Vj, and 7j(=5/3) are respectively the 
mass, density, velocity, pressure and the adiabatic poly- 
tropic index of a species j. Subscript j stands for electron 
e, proton p or alpha particle a (He 2+ ). The quantity Vj~> is 
the collision frequency of a particle of species j with par- 
ticles of species j' (only the electron-proton collisions are 
taken into account). The electric field E and the magnetic 
field B are linked by Faraday's law: 

w (6) 



V x E = - 



3.1. Linearization procedure 

The linear perturbation analysis is performed by express- 
ing all the quantities in the fluid equations as a sum of an 
unperturbed stationary part (with subscript 0) and a per- 
turbed part (with subscript 1) that is much smaller than 
the stationary part: 

rij = n 0j (z) + nij, Tj = T j(z) + Ty, pj = p 0j (z) + pij, 



■vij, B = B Q (x,z) +Bi, E = E +Ei, 



with: nij < n 0j , Ty < T Qj , pij < p j, \vij\ < |v j| , 
IBxKIBqI, |E 1 |«|E | and v 0j xB„=0. (7) 

We assume charge neutrality and a current-free state for 
the unperturbed stationary plasma, i.e., J2j1j n oj = 
and Y]j QjUoj^oj = 0. The zero-order terms cancel out 
when Eq. ^ is inserted into the multi-fluid Eqs. ([3|-([5]). 
Neglecting the nonlinear products of the first-order terms, 
we get a system of coupled linear equations: 



i(u> — k • voj)— — — ik ■ Vij = 0, 

n 0j 



(8) 



i(oj — k • v j)v 



Ei 
IBnl 



J3, 

|Bc 



■V j 



Br 

|Bo| 



v n' 



(vij 



) = o, 



■7J— = °- 
Poj n 0j 

ik x Ei = iwBi, 



(9) 

(10) 
(11) 



where all the perturbed quantities have been expressed 
in form of a plane wave. This Fourier analysis turns all 
derivatives into algebraic factors, i.e. d/dt — > — iu; and 
V — > ik., where uj is the wave frequency and k the wave 
vector. In the above equations, C S j = y/ jjkgTj /rnj is the 
acoustic speed and flj = eBo /rrij the cyclotron frequency 
of species j, and ks is the Boltzman constant. 



3.2. Dispersion relation 

To derive the dispersion relation, the above linearized 
equations have to be combined in order to obtain a linear 
relation between the current density Ji and the electric 
field Ei: 



Ji c Ei, 



(12) 



where cr is the conductivity tensor which is related to the 
dielectric tensor e through the following relation: 



e(w,k,r) = H er(cj,k,r). 

u>e Q 



(13) 



Finally, the local dispersion relation is obtained using the 
theory of electrodynamics (see, e.g., Stix 1992): 



D(w,k,r) = Dct 



-k x (k x E) +e(w,k,r) • E 



-0,(14) 



where c is the speed of light in vacuum and r is the large- 
scale position vector. We choose the wave vector k to lie 
in the x — z plane, with k = k{sin9, 0, cosO). 

3.3. Ray-tracing equations 

In the framework of the WKB approximation, the ray 
tracing problem consists in solving a system of ordinary 



differential equations of the Hamiltonian form (Weinberg 
1962). The ray-tracing equations, which represent the 



equations of motion for the wave frequency cu, the wave 
vector k, and the space coordinate r, have been formu- 



lated by Bernstein & Friedland (19841. In the simple case 



of a hcrmitian dielectric tensor, they are given by: 



duj 


dD{w,k,r)/dt 


~dt = 


dD(uj,k,r)/du 


dk 


dD(uj,k,r)/dr 


dt ~ 


&D(w,k,r)/<9w' 


dr 


&D(w,k,r)/<9k 


dT _ 


dD(uj,k,r)/duj 



(15) 



(16) 



(17) 



A generalization of these equations to the case of an anti- 
hcrmitian dielectric tensor was also proposed in the paper 



of Bernstein & Friedland (1984 1. In this case, in addition 



to the ray path, the growth rate of the instability can be 



computed as well. Note that Eq. ( 15 1 can be set to zero, 
because the dispersion relation does not explicitly depend 
on the time t (the background plasma is stationary) . The 
above set of differential equations represents an initial- 
value problem which can be solved by using the initial 
conditions obtained from the local solutions of the disper- 
sion relation (Eq. |l4"|) . 
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x = 7.5 Mm, z = 2.2 Mm, (p « 82 c 




Fig. 2. Two- fluid (e-p plasma on the left panels) and three-fluid (e-p-He 2+ plasma on the right panels) dispersion 
surfaces (top panels) and single curves (for 6 = 30° at the bottom) for the case of a zero beam speed, i.e. vo a = 0. 
The alpha-particle (Hc 2+ ) density is nQ a — 0.1no p and the plasma beta values are e — 0.0043 = 1.2/3 p = 12/3 a . Here 
u> and k are normalized, respectively, to the proton cyclotron frequency, il p , and the inertial length, Q p /Va v , where 
V J 4p=Bo / \j fiQnQpm p is the proton Alfven speed (/io is the magnetic permeability in the vacuum). Here Toe = Tq p = 



Tq 



1.35 x 10 4 if, 



nop 



3.73 x 10 16 r 



B = 15.54 Gauss. 



4. Numerical Results 

We assume that the ion-beam particles are generated at 
the funnel location at x = 7.5 Mm and z = 2.2 Mm, pre- 
sumably by small scales reconnection events. This loca- 
tion is characterized by a magnetic field inclination an- 
gle of (p = 82° with respect to the normal on the solar 
surface. According to the observations, the calculations 
will be performed for a beam velocity equal to 320 km/s. 
We consider the case of an alpha-particle (He 2+ ) beam 
plasma configuration, propagating parallel to the ambient 
field. For comparison purpose, the dispersion diagrams in 
the case of a plasma without the beam are also presented. 
We first present the results obtained from the local solu- 



tions of the dispersion relation ( 14 1 and then the results 



obtained from the non-local wave analysis using the ray- 
tracing equations. 

4.1. Local stability analysis 
4.1.1. Zero beam plasma 

For vo a — and 9 — 30°, the dispersion diagram in the 
case of the two-fluid (c-p) model (left panels of Fig. [2]) 
shows the presence of three stable modes, each one of 
them is represented by an oppositely propagating (ui > 
and u < 0) pair of waves. These modes represent the ex- 



tensions of the usual Slow (red line), Alfven (green line) 
and Fast (blue line) MHD modes into the high-frequency 
domain around u> = f2 p (= eBo/rrip), where the waves 
are dispersive (for details see Mecheri fc Marsch| 2006 ) . In 
the case of the three fluid (e-p-He 2+ ) model (right panels 
of Fig. [2]) five stable modes are present (each one repre- 
sented by an oppositely propagating pair of waves) . These 
modes are: two slow modes (yellow and red lines) , two ion- 
cyclotron modes (gray and green lines) and one fast mode 
(blue line). The waves are subject to mode conversion or 
coupling, a phenomenon associated with the appearance 
of a cut-off frequency concerning the fast mode (for details 
see |Mecheri fc Marsch|[2006 1 . 



4.1.2. Alpha particles beam 

The alpha particle beam configuration consists of three 
species: electrons (with density noe), protons (with den- 
sity riQp) and a tenuous beam of alpha particles, He 2+ 
indicated by a, with a velocity VQ a parallel to the ambi- 
ent magnetic field Bo and a density n,Q a = 0.1no p . The 
protons are considered to be at rest and the electrons are 
in motion with a velocity v$ e . Alphas and electrons sat- 
isfy thus the zero-current condition, vq c — 2(no a /no e )vo a . 
The alpha-particle-beam relative density is taken from 
in-situ observations made by the Helios spacecraft, with 
n 0a /n e ~ n 0a /n 0p = 0.05 - 0.2 (|Tu et al.|[2004|. This 
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x=7.5 Mm, z=2.2 Mm, ip m 82 c 
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kV Ap /n p 

Fig. 3. Wave frequency (black) and growth rate (red) 
in a alpha particle beam plasma (no a /no p = 0.1 and 
v 0a = 320 km/s w 1.8V]ip) versus wave vector for two 
propagation angles 6*. Note, that for = 15° (i.e., quasi- 
perpendicular propagation, top panel), the left-hand reso- 
nant ion-cyclotron mode (u>l) is excited, and that for 9 — 
65° (i.e., quasi-parallel propagation, bottom panel) the 
right-hand resonant fast mode (wr) is excited, whereby 



both satisfy the condition ( 19 1 for the anomalous Doppler 
effect. Here T 0e 



T, 



Op 



Toa- 



choice is justified since it has been argued by Fcldman 
et al. ( 1996 1 that the alpha beams might originate from 



reconnection events at the base of the expanding solar 
corona. In this case (i.e. vob ^0) Fig.[3]shows that the dis- 
persion curves are strongly modified by the appearance of 
regions of instability associated with two cyclotron beam 
modes, which in the cold-plasma case are characterized by 
the dispersion relations (Cap 1978): 



kvQ a cos(tp — 9) + f2 Q cyclotron beam 1 
kv 0a cos(tp — 9) — fl a cyclotron beam 2 



(18) 



Therefore, the value of electron-proton collision fre- 
quencies in this region of the solar atmosphere (i.e. 
v ep /Q p ssll.6 and v pe /£l p «0.006) are not large enough 
to prevent from the instabilities to develop. What colli- 
sions do is to make the electron and proton bulk velocity 
equal, but the alpha differential motion (which provide the 
free energy) remain unaffected because of the high beam 
velocity considered. 

For 9 = 15° (top panel), the instability results from 
the intersection between the ion-cyclotron mode {u>l) and 
the cyclotron beam mode 2. This region of instability is 
indicated by (L) in reference to the left-hand polarization 
characterizing the ion-cyclotron mode. Indeed, these two 



distinct and initially stable modes merge, within a certain 
range (corresponding to the red curve) of the wave num- 
ber fc, into one single unstable mode which satisfies the 
resonance condition: 



ui w kvo a cos(ip — 9) — Q, c 



(19) 



This condition corresponds to a left-hand resonant cy- 
clotron excitation of the ion-cyclotron mode through the 
anomalous Doppler effect (see, Gary 1993). As shown 



in the left panel of Fig. |4j this instability extends from 
smaller k, at propagation angles around 9 « 60° with 
small growth rate, i.e. 7 w 0.03f2 p , to higher k and 
to smaller angles of propagation with a larger growth 
rate, i.e. 7 w 0.12f2 p . Since the location at x=7.5 Mm 
and z=2.2 Mm is characterized by a Bo-inclination angle 
if = 82°, we can therefore say that as k increases this in- 
stability tends to appear at increasingly oblique propaga- 
tion angles with respect to the ambient magnetic field, and 
its growth rate tends to maximize for perpendicular prop- 
agation. On the middle panel of Fig. [4] the left-hand reso- 
nant instability is shown as a function of 9 and no a , for the 
case of Vo a = 320 km/s w 1.8Va p and a normalized wave 
number kVAp/Qp = 0.5. It is clearly seen that the growth 
rate of the instability increases with increasing no a and 
gradually is also covering a wider range of 9, but it stays 
below approximately 9 w 60° . The maximum growth rate 
is 7 « O.lfip for ti 0q = 0.2n Op and 9 ss 25°. This instability 
is also presented on the right panel of Fig. [4] as a function 
of 9 and vo a and for no a /riQ P = 0.1 and kVA p /Clp = 0.5. 
It is clearly seen that the instability has a threshold in 
the beam velocity VQ a below which it does not occur. This 
threshold depends on 9 and increases from vo a w 1.5V Ap 
at 9 ss 60° to v 0a 2V Ap for quasi-perpendicular propa- 
gation 9 w 0° (knowing that the inclination angle of Bo 
is ip w 82°). 

On the other hand, for vo a 7^ and an angle of prop- 
agation 9 = 65° (bottom of Fig. [3J, the results show the 
disappearance of the left-hand resonant instability involv- 
ing the ion-cyclotron mode and the appearance of another 
kind of instability involving the right-handed polarized 
fast mode, from which the name right-hand resonant in- 
stability is derived. This instability is indicated by (R) 
and results from the intersection of the fast mode (u>r) 
with the cyclotron beam mode 2 satisfying the resonance 
condition |l9| ) . As shown on the left panel of Fig. [5] this 
instability is mainly centered around a normalized wave 
number kVAp/^p ~ 0.5, and has maximum growth rate 
7 w 0.13f2 p at a large angle of propagation 9 w 90°. The 
right-hand resonant instability vanishes for highly oblique 
propagation with respect to the ambient field (knowing 
that the inclination of Bo in this location is ip » 82°), and 
the wave amplitude grows strongly for decreasing obliquity 
of the propagation angle. On the middle panel of Fig. [5] 
the right-hand resonant instability is shown as a function 
of 9 and n 0a , for the case of v 0a = 320 km/s w 1.8Va p 
and kVAp/Qp — 0.5. It is clearly seen that the growth 
rate of this instability increases with increasing no a - This 
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x = 7.5 Mm, z = 2.2 Mm, ip sa 82° 



voa/Vip 1.8, noa/nop = 0.1 v 0a /V Ap ss 1.8, feF^/Oj, = 0.5 ra Oa /«0j> = 0.1, feVip/fij, = 0.5 




kV Ap /n p QoJ%p Vo^/Vap 

Fig. 4. Growth rate of the left-hand resonant ion-cyclotron instability, in the case of an alpha-particle beam 
plasma. Left: as a function of the angle of propagation 9 and the normalized wave number kVA p /Q p and for an alpha 
beam density no a — 0.1no p and velocity VQ a — 320 km/s = 1.8Va p - Middle: as a function of 9 and no a and for 
VQ a = 320 km/s w 1.8V^ P and fcVA p /fi p = 0.5. Right: as a function of and voq/wa p and for riQ a — 0.1no p and 
kV Ap /£l p = 0.5. Here T 0e = T 0p = T 0a . 



behavior is more pronounced for higher 9, which corre- 
sponds to decreasingly oblique propagation. The maxi- 
mum growth rate, 7 w 0.16f2 p , is obtained for no a =0.2no p 
and 9 « 80° . We can also notice that the instability fades 
away for small propagation angles, with 9 < 40°, which 
corresponds to quasi-perpendicular propagation (with re- 
spect to the field) . In the same figure (on the right panel) 
we also present the dependence of this instability upon 9 
and v Qa for no a /no p = 0.1 and kV Ap /£lp — 0.5. It can 
be seen, similarly to the left-hand instability, that this in- 
stability has a threshold in the beam velocity vo a , below 
which it docs not occur. This threshold depends on 9 and 
increases from VQ a rs 1.6V a p at large angle of propaga- 
tion, 9 f=s 90° (quasi-parallel propagation since ip — 82°), 
to Vo a « 2.5Vyi p at 9 ~ 25° (quasi-perpendicular propa- 
gation) . 

4.2. Non-local stability analysis 

In this section we intend to go beyond the local treatment 
of the waves and perform a non-local wave study using 
the ray-tracing equations. The ray-tracing equations are 
solved employing the initial conditions obtained from the 



local solutions of the dispersion relation ( 14 1 at the loca- 
tion with xq = 7.5 Mm and zq = 2.2 Mm. We consider 
the case of an alpha-particle (He 2+ ) beam plasma config- 
uration with a constant concentration, n a = 0.1n p , and a 
constant beam velocity of vo a = 320 km/s. 

The ray paths of the unstable waves as well as the vari- 
ation of their growth rates as a function of height z, when 



the wave is launched at the initial location (xq — 7.5 Mm 
and zq = 2.2 Mm) that is characterized by a strong in- 
clination angle (ip as 82°) of the magnetic field with re- 
spect to the normal on the solar surface, are illustrated in 
Fig. [6j The results are presented for a different initial an- 
gle of propagation, #0, with which an initial wave number 
ko (that is normalized to £I p /Va p ) is associated and chosen 
as to correspond to the maximum growth rate, 7 mQX . 

Our results show that the ray path of the left-hand 
unstable wave (Fig. [6] on the left) is strongly affected by 
the closed-field geometry characterizing this funnel region. 
Indeed, this unstable wave starting from its initial posi- 
tion propagates upward in the coronal funnel to a certain 
height, where it turns down again and starts propagating 
downward to return back to the initial height. The as- 
sociated instability growth rate decreases along that ray 
path. Since the direction of the group velocity is always 
parallel to the ray path and indicates where the energy is 
transported, we can say that the energy associated with 
the left-hand resonant instability does not reach high al- 
titudes in the funnel. The smaller 9q is the higher up this 
unstable wave propagates in the funnel. 

The right-hand unstable wave (Fig. [6] on the right) is 
also found to be well guided and, depending on #0, can 
follow both closed and open coronal field lines, which may 
exist side by side in this region of the funnel. Indeed, for 
9 = 50° and #0 = 55° the unstable wave propagates along 
the open magnetic field lines and reaches high altitudes in 
the funnel up to 15 Mm, while for 9q = 60°, 65° and 
70°, similarly to the left-hand instability, the right-hand 
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v 0a /V Ap m l,8, n 0a I n 0p = A »oi*/Va s ~ 1.8, kV^ p /r! p =0.5 n 0a jn 0v = 0.1, kV 4 j,/fij,=0.5 




Fig. 5. Growth rate of the right-hand resonant ion-cyclotron instability, in the case of an alpha-particle beam 
plasma. Left: as a function of the angle of propagation 9 and the normalized wave number kVAp/Qp and for an alpha 
beam density no a — 0.1no p and velocity vo a — 320 km/s = 1.8Vyip. Middle: as a function of 9 and no a and for 
voa = 320 km/s w 1.8Vyip and kVAp/Qp = 0.5. Right: as a function of 9 and vob/vA P and for no Q = 0.1no p and 
kVAp/Qp = 0.5. Here T 0e = T 0p = T 0a . 



resonant unstable waves are affected by the closed-field 
geometry and reflected back towards lower altitudes in 
the funnel. Thus, the energy associated with this insta- 
bility is also transported along the magnetic field lines, 
but eventually to much greater altitude, i.e. z = 15 Mm, 
in the funnel as compared to the left-hand instabilities. 
For 9q = 50° and 9q = 55° the growth rate of the right- 
hand instability first increases until respectively altitudes 
of z ps 3.25 Mm and z w 3.75 Mm where it starts to 
sharply decrease and cancels respectively at z ss 3.6 Mm 
and z « 4.5 Mm. On the other hand, for 9 = 60°, 65° 
and 70°, the growth is first slightly increases during the 
upward propagation phase, and then it rapidly decreases 
to a zero value while the wave is propagating downward. 



5. Conclusion 

We have studied beam-driven electromagnetic instabili- 
ties near the ion-cyclotron frequency in a coronal fun- 
nel using the multi-fluid model. We have considered the 
case of an alpha-particle beam propagating in the fun- 
nel parallel to the ambient magnetic field lines. In agree- 
ment with kinetic dispersion theory, the local solutions 
of the dispersion relation revealed the presence of two 
kinds of instabilities: the left-hand and right-hand reso- 
nant instabilities. The left-hand and right-hand instabili- 
ties arise from the resonant excitation of, respectively, the 
left-hand-polarized ion-cyclotron mode and right-hand- 
polarized fast mode through the anomalous Doppler ef- 



fect, see Eq. (19). For the studied coronal region, our re- 



sults indicate that the left-hand resonant instability devel- 



ops for strongly oblique wave propagation with respect to 
the ambient magnetic field, with a maximum growth rate 
at a quasi-perpendicular propagation, and disappears for 
weakly oblique propagation. Oppositely, the right-hand in- 
stability develops for a weakly-oblique propagation to the 
ambient magnetic field, with a maximum growth rate at 
quasi-parallel propagation, and ceases for highly oblique 
or quasi-perpendicular propagation. 



The nonlocal ray-tracing analysis revealed that both 
instabilities are sensitively affected by the magnetic field 
geometry and found to propagate closely along the field 
lines. The left-hand resonant instability is rapidly re- 
flected, thus obeying the constraints imposed by a closed- 
field configuration. The associated growth rate slightly de- 
creases and eventually cancels along the ray path. On the 
other hand, the right-handed resonant instabilities are also 
very well guided along the magnetic field lines of the fun- 
nel. This instability, for small initial angles of propaga- 
tion, appears to follow the open field lines and can propa- 
gate higher up in the funnel, yet with a rapidly decreasing 
growth rate. 



Consequently, fast ion beams in the magnetically open 
corona can provide enough energy for driving micro- 
instabilities through resonant wave-particle interactions. 
These instabilities may constitute in turn an important en- 
ergy source for high-frequency ion-cyclotron waves which 
have been invoked to play a relevant role in the heating of 
coronal ions through cyclotron damping. 
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Fig. 6. The ray trajectory (top) and growth rate (bottom) of the left-hand (left) and the right-hand (right) resonant 
ion-cyclotron instabilities in the coronal funnel. These unstable modes are launched at the initial location (xo = 7.5 Mm 
and zq — 2.2 Mm) in the funnel (with a Bo-inclination angle <po f=a 82°), for different initial angle of propagation 0q 
and wave number k chosen at the maximum instability growth rate. The dashed lines represent the magnetic field 
lines in the funnel. 
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